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ABSTRACT 



Context. Extragalactic jets are formed close to supermassive black-holes in the center of galaxies. Large amounts of gas, 
dust, and stars cluster in the galaxy nucleus, and interactions between this ambient material and the jet base should 
be frequent, having dynamical as well as radiative consequences. 

Aims. This work studies the dynamical interaction of an obstacle, a clump of matter or the atmosphere of an evolved 
star, with the innermost region of an extragalactic jet. Jet mass- loading and the high-energy outcome of this interaction 
are briefly discussed. 

Methods. Relativistic hydrodynamical simulations with axial symmetry have been carried out for homogeneous and 
inhomogeneous obstacles inside a relativistic jet. These obstacles may represent a medium inhomogeneity or the dis- 
rupted atmosphere of a red giant star. 

Results. Once inside the jet, an homogeneous obstacle expands and gets disrupted after few dynamical timescales, 
whereas in the inhomogeneous case, a solid core can smoothen the process, with the obstacle mass-loss dominated by a 
dense and narrow tail pointing in the direction of the jet. In either case, matter is expected to accelerate and eventually 
get incorporated to the jet. Particles can be accelerated in the interaction region, and produce variable gamma-rays in 
the ambient matter, magnetic and photon fields. 

Conclusions. The presence of matter clumps or red giants into the base of an extragalactic jet likely implies significant 
jet mass-loading and slowing down. Fast flare-like gamma-ray events, and some level of persistent emission, are expected 
due to these interactions. 

Key words. Galaxies: jets-Gamma rays: galaxies-Stars: AGB and post-AGB 
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1. Introduction 

Jets of active galactic nuclei (AGN) are collimated outflows 
originated at the core of galaxies, very likely in the vicinity 
of a supermassive black hole (SMB H) that accret es matter 
from its environment (e.g. Beg elman et all Il984l ). The in- 
ner regions of galaxies contain large amounts of gas, dust, 
and stars (e.g. iBurbidgel [l970), so there is plenty of ma- 
terial that could interact with the AGN jet. About a 1% 
of th e present stars are evolved objects (e.g. I Young et all 
Il977f ) whose external layer s, close to the SMB H, can be 
tidally distorted (see, e.g., iBarkov et all I2Q1QL and refer- 
ences therein). Stellar c ollisions coul d also release matter 
to the environment (e.g. lYoungl Il977l ) . In AGN with signif- 
icant star formation in their cores, massive stars with strong 
winds will be also present. In bright AGN, hot gas embed- 
ding denser and cooler clouds is expected to surround the 
SMB H (th e Bro ad Line Region -BLR-; e.g. iKrolik et all 
ll98ll : lReeslll987[ ). In these bright AGN, X-ray heating can 
also perturb the e xterna l layers of stars even i f tidal forces 
are nee; ligible (e.g. iShullL fl98l iPenstonl . fl988). Wandering 
dark clouds could also be present in the region close to 
the SM BH, as it has been re cently discovered in our own 
galaxy (|Gillessen et all l2Q12h . It seems therefore unavoid- 
able that stellar matter and medium inhomogeneities will 



frequently interact with, and in some cases penetrate into, 
the AGN innermost jet regions. Entrainement of medium 
matter leads to jet mass- loading, invoked by seve ral authors 
(e.g. lBowman et all 11996c lLaing fc Bridlel 12002) to explain 
the deceleration of FRI jets between pc-scales, where they 
show relativistic speeds, and kpc-scales, where their veloci- 
ties are sub-relativistic as shown by the jet and counter-jet 
symmetric emission. 

Studies of the dynamics of the interaction between 
clouds or red giants (RG) after e ntering into an AGN 
jet have been carried out (e.g. lAraudo et all l2010b 
I Barkov et alll2Q10Ll2QTTh , However, previous works were of 
analytical nature, mainly focused on the subsequent high- 
energy emission. Other works have studied the global im- 
pact of stars or clouds on the jet propagation and con- 
tent ( e.g. iKomissarovl 1 19941: ISteffen et all 11997c IChoi et al ' 



20051: Hubbard fc BlackmanL 120061 : ISutherland et all 12007 



Jevakumarl l2QQ9h . but numerical studies of the conse 
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quences of obstacles being inside the innermost regions of 
AGN jets are rare. 

In this work, we aim at studying numerically the evo- 
lution of a dense cloud and the disrupted layers of an RG 
inside an AGN jet. Homogeneous and power-law (plus a 
solid core) density profiles have been adopted, the former 
to simulate the simplest cloud or RG scenario, and the lat- 
ter what might be a more realistic situation for a disrupted 
RG atmosphere. This can allow a characterization of the 
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dynamical evolution of obstacles inside jets more accurate 
than in previous works, as well as a better understand- 
ing of the impact on the jet itself, and the resulting high- 
energy radiation. At this stage, we neglect the role of the 
magnetic field in the jet, accounting only for its ram pres- 
sure, and assume an axisymmetric (2D) interaction, with 
the star/obstacle at rest in the laboratory frame. Magnetic 
fields and 3-dimensional (3D) effects will be included in fu- 
ture work. 



2. Physical scenario 

Unlike smooth media, which interacts with the jet through 
a shear layer preventing direct mixing, medium inhomo- 
geneities like matter clumps, stars, etc., could penetrate 
effectively into the jet. The requirement for the obstacle 
to enter into the jet is to have a jet perpendicular speed 

^ ^sc, where v sc is roughly the shocked obstacle sound 
speed, i.e. v sc ~ Cy / pfi-J p Q , T } = 1/(1 - (v } /c) 2 ) 1/2 the jet 
Lorentz factor, and pj and p Q the jet and obstacle densities 
in the reference frame of the latter. This assumes that the 
relevant jet pressure comes from bulk motion, which im- 
plies that the jet power Lj = 7rR 2 (Tj — l)pj^jC 2 , where R } is 
the jet radius. Once the obstacle is inside the jet, neglect- 
ing magnetic field effects, a bow-shaped shock forms in the 
jet. Another shock with speed v sc crosses t he obstacle, or 
in th e case of an RG , its ex ternal layers (e.g. lAraudo et aH 
201(J; IBarkov et al.L I2Q1QL and references therein). Being 
much faster, the jet is to be much lighter than the ob- 
stacle to entrain the latter. Therefore, the obstacle shock 
crossing (or dynamical) time td ~ R /v sc , with R Q as the 
obstacle radius, will be much longer than the bow-shock 
formation time, ~ R /c. After few td, the obstacle has 
been accelerated by the jet up to a vel o city ~ v\ (e.g. 
iBlandford fc KoeniglL 119791 : IBarkov et all I2010L l2011h . It 
follows that if the time of jet crossing, tj ~ R]/v Ql is smaller 
than td, the obstacle will be able to leave the jet; otherwise, 
it will be dragged by the jet and presumably disrupted and 
mixed with its matter. 

For obstacles remaining long enough in the jet, and 
to have a dynamically strong impact on the jet, the ob- 
stacles have to affect significantly the total jet mass, mo- 
mentum and energy fluxes. In the present context, it can 
mean also that enough mass will be entrained to enforce 
significant jet deceleration and energy dissipation, i.e. the 
obstacle mass-injection flux M ~ M cr = Z/j /(Tj — 1) c 2 
([Hubbard fc Blackmanl l2QQ6f ). In addition, tj should be 
significantly longer than ^, to allow the jet to drag the 
whole obstacle matter. The value of M can be estimated 
as ~ A/q x 2 M /4tj(z ), where N Q is the number of ob- 
stacles (clouds, stars, etc.) within a sphere of radius z Q , 
the jet characteristic height, x — R}( z o)/zo °° 0-1, and 
M Q = 47ri?oPo/3 the detachable obstacle mass. The value 
of z Q depends strongly on th e scenario under stu dy. It 
might be the size of the BLR (Ar audo et aH . I201QD . a re- 
gion in which SMBH tidal forces affect the RG atmosphere 
([Barkov et all |201Q[ ), or the distance at which accretion 
disk X-rays can affect the star external layers. Note that 
for td <C tj, the obstacles will initially affect just the outer- 
most jet shell, but eventually a strong shear layer will de- 
velop reaching the jet core. For td ^ tj, the obstacle mass 
will be more distributed over the jet cross section, likely 
speeding up the process of mass-loading and energy dissi- 



pation (e.g. emission). Note that here only sudden mass- 
loss by shocked obstacles is treated; continu ous mass-loss 
via, e.g., stellar winds is not considered (see I Araudo et all 
120121 ; Perucho et al., in preparation). A significant dynam- 
ical impact will likely have its observable counterpart in 
the form of non-thermal radiation. However, we remark 
also that, even when providing little mass- loading, the oc- 
casional penetration of dense and large o bstacles within 
the jet could have observable co nsequences ([Araudo et all 
120101 : IBarkov et all I2010L l2011h . 

Figure [1] shows a sketch of the scenarios simulated here. 
The simulations account for the phase when the obstacle is 
well inside the innermost region of an AGN jet. Two types 
of obstacles have been considered. The first type consists on 
an homogeneous gas sphere, which may represent a dense 
cloud in the central AGN regions (e.g. from the BLR), or a 
homogeneous bulk of material from the external RG layers 
detached by tidal SMBH forces. The second case considers 
a power-law profile with a radial dependence ex R~ 2 for the 
density (plus a solid core), which may represent the RG 
external layers more realistically, or just slightly affected 
by tidal forced]. 




Accretion disc 



Fig. 1. Sketch of the two scenarios contemplated here, a) 
Penetration of an RG, with the external layers detached (to 
different degrees) due to gravitational disruption, into the 
jet. b) Penetration of a massive clump of matter into the 
jet. 



3. Simulations 

We have performed three numerical simulations. In simu- 
lation 1 (SI), we have calculated the interaction between 
a homogeneous spherical obstacle and the jet; in simula- 
tion 2 (S2) we have introduced a constant density core plus 
a power-law profile for the radial variation of the density 
out of this core (R~ 2 , see Section [2]); in simulation 3 (S3) 
we have reproduced S2 with double resolution. The three 
simulations have 2D axial symmetry and have been carried 
out in cylindrical coordinates. This is justified by the sym- 
metric nature of the problem solved (initial interaction and 

1 It is worth noting that the power-law case may also cor- 
respond to an RG with its wind confined by external pressure 
outside the jet, forming a bubble that would be brought by the 
RG into the jet. 
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obstacle ablation), which may only be broken by inhomo- 
geneities in the jet that lead, at larger scales than studied 
here, to instabilities in the cometary tail formed by the ob- 
stacle material and subsequent mixing with jet matter. It 
should be noted that these simulations represent a simplifi- 
cation, and a necessary initial step, of the more realistic sit- 
uation, which would require the obstacle to penetrate into 
the jet with transversal motion. This can only be solved via 
3D-simulations and is beyond the scope of this initial work. 

We have used the finite-volume code Ratpenat, which 
solves the equations of relativistic hydrodynamics in con- 
servation form using high-resolution-shock-capturing meth- 
ods. Ratpenat is a hybrid parallel code - M PI + OpenMP 
- ext ensively and intensively tested (e.g., iPerucho et al.L 
201(j). The c ode includes the Synge equation of state 
(jSvngel fl957h with two populations of particles, namely, 
leptons (electrons and positrons) and baryons (protons). 
The simulations were performed in Tirant, at the Servei 
dinformatica of the Universitat de Valencia, using 16 pro- 
cessors in SI and S2, and 32 processors in S3. 

In SI, the numerical grid box expands to 10 15 cm in the 
radial direction, with a homogeneous grid up to 2 x 10 14 cm 
including 400 cells, plus an extended grid with 200 cells 
covering the remaining 8 x 10 14 cm. In the direction of the 
jet flow (axial direction), the grid covers 6 x 10 14 cm. The 
(homogeneous) obstacle is located at z = 10 14 cm, with 
radius R Q = 10 13 cm, and the rest of the grid is filled by 
the jet flow. The resolution used in SI is 20 cells/it^ with 
a total of 600 x 1200 cells. In S2, the grid expands to 3 x 
10 14 cm, with a homogeneous region up to 10 14 cm using 320 
cells, and 200 cells extending the grid 2 x 10 14 cm further. 
Along the symmetry axis, the grid is 3 x 10 14 cm long. The 
obstacle is located at z = 2.5 x 10 13 cm, it has a radius R Q = 
10 13 cm, with an inner core of constant density with radius 
of 2.5 x 10 12 cm, and a surrounding region with decreasing 
density (see Sect. [2|) up to R Q . As in SI, the rest of the 
grid is filled by jet flow. The resolution used in S2 is 32 
cells/i^o, with a total of 520 x 960 cells. In the case of S3, 
the simulation is focused on a smaller region of S2, using a 
resolution 64 cells/it^ with a total of 320 x 960 cells that 
cover a grid size of 5 x 10 13 cm in the radial direction and 

1.5 x 10 14 cm in the axial direction. 

The jet is modelled in all three simulations as a mat- 
ter dominated flow with proton-electron composition. The 
properties of the jet in SI have been chosen to result in 
a total jet kinetic power Lj >s i = 2 x 10 44 erg/s for a jet 
radius Rj = 10 15 cm, with jet velocity v } = 0.866 c = 

2.6 x 10 10 cm/s, density p- hsl = 1.35 x 10" 18 g/cm 3 , tem- 
perature Tj = 10 10 K, adiabatic exponent 7j = 1.45, and 
Mach number Mj = 16.8. With these parameters, the ki- 
netic power of the jet within the homogeneous simulated 
grid is £j, g i = 8 x 10 42 erg/s and the mass flux is Mj, g i = 
2.8 x 10 21 g/s. In S2 and S3, the jet density used is slightly 
larger />j, S 23 = 2.4 x 10 _18 g/cm 3 , giving a total jet power 
within a radius Rj of £j, S 23 = 3.5 x 10 44 erg/s. In the case of 
S2, Lj, g2 = 3.5 x 10 42 erg/s and M hg2 = 1-25 x 10 21 g/s, and, 
for S3, Lj, g3 = 8.75 x 10 41 erg/s and M- hg3 = 3.12 x 10 20 g/s. 

The obstacle in SI has density p Q = 1.35 x 10 _12 g/cm 3 
and temperature T Q = 10 4 K, with a total mass M Q ~6 x 
10 27 g. In the case of S2, the obstacle has a core density 
p = 1.25 x 10 _10 g/cm 3 within the inner 2.5 x 10 12 cm, 
and a decreasing density ex R~ 2 (R is the spherical radius 
here), starting with p = 1.7 x 10 -11 g/cm 3 at 2.5 x 10 12 cm. 



The material of the shocked obstacle can have rather 
complex properties. Once in pressure equilibrium with the 
jet ram pressure, the obstacle is radiation dominated, and 
radiative cooling through thermal Bremsstrahlung may be 
important. For simplicity, we have modelled the obstacle 
as a non-relativistic proton-electron plasma. This has some 
implications for the compression and stability of the ob- 
stacle, since the adiabatic coefficient may not be accurate 
with respect to realistic cases, and shocked matter densities 
(and instabilities) could be enhanced by radiative cooling. 
However, at this stage of the study the calculations can 
capture the relevant aspects of the simulated process. More 
accurate simulations will be performed in the future. 

4. Results 

4.1. Homogeneous case 

In SI, the interaction region rapidly develops the expected 
structure, with a bow shock in the jet, and a contact discon- 
tinuity separating both the shocked media and the obstacle 
shocked by a slow shock. The eroded material from the 
obstacle is dragged downstream, forming a cylindrical tail 
surrounded by rarefied jet flow that has crossed the bow 
shock. 

Figure [2] shows a snapshot of the simulation at t c± 9 x 

10 4 s, which corresponds to the initial phase. After t = 2.3 x 

10 5 s, the shock has completely crossed the obstacle. The 
pressure at the apex of the bow shock oscillates between 10 3 
and 1.5 x 10 3 dyn/cm 2 , and the temperature of the post- 
shock gas is at this point ~ 10 12 Kq After the obstacle 
matter is heated by the shock, it expands, increasing the 
cross-section of the interaction with the jet flow. Up to this 
point, the mass- load by obstacle material is modest (see 
Sect. 14. 3p and mainly due to continuous ablation from its 
outer layers. This gas accumulates in a cylindrical region 
with radius of the order of R , forming a cometary-like tail 
of obstacle gas separated from the rarefied, shocked jet flow 
by a smooth shear layer in which little mixing occurs, as no 
instabilities are triggered in the layer during this phase. 

Figure [3] shows the trends in energy flux, mass flux, ob- 
stacle mean density (averaged over the grid area Aj) and 
velocity along the axial direction for a snapshot at a rela- 
tively early stage. In this plot, the energy flux is separated 
into kinetic plus rest-mass energy flux (£j, g i — £j,int,gi = 
(p(hW — 1)W - pWs)vjAj, with h being the enthalpy, W 
the Lorentz factor, and e the specific internal energy), and 
internal energy flux (Lj,i n t,gi = pWe)v- } A- } ). These variables 
give an idea of the amount of energy flux in the form of in- 
ternal energy or kinetic and rest-mass energy. A fraction of 
Lj )g i in this phase is mainly invested in heating the shocked 
jet material, with only a small amount transferred to the 
obstacle gas (upper panels in Fig. [3j). Some fraction of this 
energy may go to non- thermal particles. The gas in the tail 
expands and accelerates along the z coordinate, reaching a 
mean axial speed v z ~ 1.2 x 10 8 cm/s (bottom left panel in 
Fig. [3]), with v z ~ 10 9 cm/s right on the symmetry axis. 

After the obstacle is shocked, instabilities develop inside 
and disrupt it, detaching clumps of material and enhancing 
turbulent mixing with the jet flow in the cometary-like tail. 
The instabilities manage to tear the external portions of ob- 
stacle gas at t ~ 7 x 10 5 s, as the jet flow penetrates through 

2 The tem perature could be significantly lower due to pair 
production (e.g. iBisnovatvi-Kogan et allll97lh . 
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Fig. 2. Combined maps of logarithm of pressure (left panel, upper half), density (left panel, lower half), axial velocity 
(right panel, upper half), and jet mass fraction (right panel, lower half) at t c± 9 x 10 4 s for SI. The plots make use of 
the axisymmetric nature of the simulation. Jet mass fraction equals 1 for pure jet material, for pure obstacle material, 
and values in between indicate mixing. Units are cgs. 
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Fig. 3. Axial cuts of different quantities for SI at t ~ 6 x 10 5 s, before the instabilities completely destroy the obstacle. 
The top left panel shows Lj jg i — £j,int,gi for the jet material (solid line) normalized to its injection value 8 x 10 42 erg/s, 
and the internal energy luminosity Lj,i n t,gi normalized to its maximum value, 8 x 10 40 erg/s (dashed line); the red solid 
line stands for the internal energy luminosity in the jet material, and the blue dotted line (very close to the bottom) 
for the internal energy luminosity in the obstacle material. The top right panel shows the mass flux, normalized to its 
maximum value (4.7 x 10 22 g/s). The bottom left and right panels show the mean velocity (averaged over the grid section) 
and density in the obstacle gas. 

the obstacle. This process repeats itself in smaller clumps rial remaining at the symmetry axis, all of them also torn by 
that have been detached from the obstacle and in the mate- the jet (see Figs. [Hand [5]). These clumps move radially and 
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Fig. 4. Same as Fig [2] at t~ 1.1 x 10 6 s 



are accelerated when entering the region of rarefied (post- 
shock) jet material, reaching mean velocities ~ 3x 10 9 cm/s. 
At the same time, the jet ram-pressure starts to effectively 
push the obstacle remains from its original location along 
the z-direction, dragging them downstream to the end of 
the grid. The left and central panels of Fig. [6] show the in- 
crease in mass flux and mean density when the rest of the 
obstacle crosses through the axial locations where these val- 
ues are computed (half and end of the grid, along the axis). 
At the same time, the right panel shows a strong decrease in 
the mean velocity at the plotted locations, which is accom- 
panied by a drop of the mean Mach number from Mj = 16 
to values < 5. 

Figure [7] shows the same plots as Figure [3] but at 
t ~ 1.15 x 10 6 . At this time the cross-section of the obstacle 
is maximum, as shown by the upper panel in Figure [5] The 
plots show that the transfer of energy from the jet material 
to the obstacle is maximal in this situation, with the jet 
losing almost all of the initial kinetic and mass energy flux, 
which is put into internal energy flux of the jet and the 
obstacle (and potentially into non-thermal particles), plus 
kinetic and mass energy flux of the obstacle material (not 
shown). The latter is revealed by the large mean velocities 
reached by the obstacle material (bottom left panel). The 
mass flux shows an increase by nearly two orders of mag- 
nitude with respect to the original jet flux in the simulated 
region, and the obstacle mean density along the grid is also 
larger by several orders of magnitude, depending on the 
position along the axis, than that shown in Figure [3] 



4.2. Inhomogeneous case 

In this case, the interaction of the jet with the obstacle 
is steadier than in the homogeneous case. The pressure 
at the apex of the bow shock oscillates between 2 x 10 3 
and 2.5 x 10 3 dyn/cm 2 , and the temperature of the post- 
shock gas there is again ~ 10 12 K. The ablation of mass 
from the outer, more dilute layers of the obstacle occurs 
at a smoothly increasing rate, until t ~ 8 x 10 5 s, when 
a significant portion of those layers is detached from the 
denser core. After that, almost the whole obstacle beyond 
2.5 x 10 12 cm has been shocked, and the tail is loaded with 
this gas. At obstacle radii < 2.5 x 10 12 cm the density is 
so high that the material is not affected by the jet during 
the simulation (as expected for the stellar layers below the 
photosphere). Unlike in SI, shocking the obstacle external 



% 2.03e-16 



% 1.149-16 




Fig. 5. Same as left panel in Fig. [2] at t ~ 1.15 x 10 6 s, 
t ~ 1.3 x 10 6 s and t - 1.45 x 10 6 s. 
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Fig. 6. Mass flux (left panel), mean density (central panel), and mean velocity (right panel) versus time for SI at 
z = 3 x 10 14 (half grid, black lines) and 6 x 10 14 cm (end of the grid, red lines). 
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Fig. 7. Axial cuts of different quantities for SI at t c± 1.15 x 10 6 s, the time of maximum cross-section of the interaction 
between the obstacle and the jet (see Fig. [5j). The top left panel shows Lj >g i — £j,int,gi for the jet material (solid line) 
normalized to its injection value 8 x 10 42 erg/s, and the internal energy luminosity Lj,mt,gi normalized to its maximum 
value, 3.6 x 10 42 erg/s (dashed line); the red solid line stands for the internal energy luminosity in the jet material, and 
the blue dotted line for the internal energy luminosity in the obstacle material. The top right panel shows the mass flux, 
normalized to its maximum value (10 23 g/s). The bottom left and right panels show the mean velocity and density in the 
obstacle gas. 



layers does not lead to an enlargement of the bow shock 
cross-section, since the radial drop in density smoothens 
the process to some extent. As shown later, this effect is 
enhanced by a low resolution (see Sect. I4.2.T]) . 

Figure [8] shows a series of maps of the whole process. 
The upper left panel corresponds to the initial phase. The 
next two panels (top right and central left), show the shock 



travelling around the core and through the outer layers of 
the obstacle. In the central right panel, these layers have 
already been completely shocked and a large amount of gas 
that has been detached from the core is visible as a bump in 
the tail between z = 5 x 10 13 cm and z = 7 x 10 13 cm. The 
bottom panels show this bulk of mass propagating down- 
stream. 
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Fig. 8. Same as left panel in Fig. for simulation S2 at times t ~ 5 x 10 4 s, t ~ 1.6 x 10 5 s, t ~ 4.3 x 10 5 s, t ~ 7.9 x 10 5 s, 
t ~ 1.25 x 10 6 s and t ~ 1.8 x 10 6 s from top to bottom and left to right. 



The upper panels in Fig. [9] show the normalized Lj )g 2 — 
£j,int,g2 and mass flux versus distance along the axis at 
t ~ 3.1 x 10 5 s, i.e. still during the initial phase. As seen 
in the figure, the jet loses a small amount of kinetic energy 
flux, which is used to heat the jet material and to heat and 
accelerate obstacle material. The mass flux falls a factor 
of ten from the obstacle itself to its tail, where it has a 
constant value. The lower panels in Fig. [9] show the values 
of mean velocity and density of the obstacle gas versus z 
at the same time. The values obtained for the density in 
this initial quasi-steady regime of S2 are similar to those 
obtained for SI. Regarding velocity, the obstacle material 
is accelerated to v ~ 1.9 x 10 8 cm/s at the end of the grid, 
which is a factor of 3 larger than the velocity reached by 
the obstacle material in SI at the same axial distance. This 
difference is in part due to the smaller total mass flux in 
S2 (see left panels in Figs. l6l and fTU]) , caused by the lower 



density of the external layers. On the symmetry axis, the 
maximum axial velocity is v z ~ 10 9 cm/s at the end of the 
grid, as in the case of SI, but at z = 3 x 10 14 cm, i.e., half 
the distance needed in SI to reach the same velocity. 



Figure [TO] shows the mass flux within the grid and the 
mean density for two different z- values, namely at half (z = 
1.5 x 10 14 cm) and at the end of the grid (z = 3 x 10 14 cm), 
versus time. These plots show that the mass load grows 
smoothly with time at z = 1.5 x 10 14 cm until t ~ 1.4 x 
10 6 s, which corresponds to the time of passage of the clump 
of mass ablated around t ~ 8 x 10 5 s. By the end of the 
simulation, this clump has not still reached the end of the 
grid. The mean velocity of the simulated portion of the jet 
is not significantly reduced during the whole simulation. 
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1.2i I 1.0 




z (cm) z (cm) 

Fig. 9. Axial cuts of different quantities for simulation S2 at t ~ 3.1 x 10 5 s, during the initial phase, before the shock 
completely crosses the obstacle. The top left panel shows Lj )g i — £j,int,gi for the jet material (solid line) normalized 
to its injection value 3.5 x 10 42 erg/s, and the internal energy luminosity Lj,i n t,gi normalized to its maximum value, 
5.6 x 10 40 erg/s (dashed line); the red solid line stands for the internal energy luminosity in the jet material, and the 
blue dotted line for the internal energy luminosity in the obstacle material. The top right panel shows the mass flux, 
normalized to its maximum value (4 x 10 22 g/s). The bottom left and right panels show the mean velocity and density 
in the obstacle gas. 




Time (s) Time (s) 

Fig. 10. Mass flux (left) and mean density (right) versus time for simulation S2 at z = 1.5 x 10 14 (half grid, black lines) 
and 3 x 10 14 cm (end of the grid, red lines). 

4.2.1. Inhomogeneous case: high resolution 

Simulation S3 reproduces a region of S2 with double resolu- 
tion, and half of its grid size along and across, resulting in 
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a physical grid size of 5 x 10 13 cm x 1.5 x 10 14 cm. From the 
results, we see that the apex pressure and temperature in 
S3 are very similar to those in S2. Figure [TT1 shows different 
snapshots of the simulation, at the same instants as for S2 
but stopping at t ~ 1.25 x 10 6 s. At first sight, both simu- 
lations follow the same phases, i.e., a quasi-steady phase in 
which the outer layers of the obstacle are being crossed by 
the shock and the mass-loss is smoothly increasing, and a 
second phase, starting when the outer region has been com- 
pletely crossed by the shock, and a large amount of gas is 
detached from the core and dragged downstream. However, 
the main difference, which can only be related to the numer- 
ical resolution, is that this detached clump of gas is larger 
and produces some widening of the jet-obstacle interaction 
cross-section (compare Figs. [8] and [11] at t c± 7.9 x 10 5 s), af- 
ter the whole outer layers of the obstacle have been shocked 
at t ~ 4 x 10 5 s. Thus, this material receives a larger transfer 
of energy from the jet, and is accelerated in a more efficient 
way, reaching z = 1.5 x 10 14 cm at t ~ 10 6 s, compared to 
t ~ 1.4 x 10 6 s in S2, as shown in Fig. [B] 

Figure [12] shows the mass flux across the simulated grid, 
and the mean obstacle gas velocity and density, versus dis- 
tance along the axis at a time similar to that in Fig. [9] We 
have omitted the luminosity plot because it shows the same 
behavior in both S2 and S3. However, at the time of maxi- 
mum obstacle expansion in S3, i.e. t ~ 7.9 x 10 5 s, the repro- 
cessed luminosity is smaller than Slatt^l.l5xl0 6 s, but 
significantly larger than in S2. Regarding the mass flux, we 
observe that it keeps smaller relative (and absolute) values 
in S3 than in S2. The tail is narrower in S3 than in S2. The 
mean obstacle gas density along the tail has similar values 
in both simulations, as well as the temperature. However, 
the mean obstacle gas velocity presents larger velocities by 
a factor of 2 in S3 than in S2. The smaller value of the 
obstacle mass-loss during the initial phase could be respon- 
sible for this difference, as indicated by the mean velocity 
of the obstacle gas right at the region of interaction, which 
is a factor two larger in S3 than in S2. 

Figure [13] shows the mass flux, mean density and mean 
velocity at z = 7.5 x 10 13 and 1.5 x 10 14 cm versus time. 
Note that a direct comparison between S2 and S3 here can- 
not be done, as the physical cross section of the numerical 
grid is different in both simulations. However, Fig. [13] al- 
ready gives interesting information regarding the mass-load 
process. In S3, the mass flux, which should be similar to 
that in S2 when dominated by the outer material of the 
obstacle, presents smaller values during the initial quasi- 
steady phase and larger values when the aforementioned 
large clump is detached from the core. This implies that the 
obstacle mass-loss is slower in S3 during the initial phase 
and, thus, after the obstacle has been completely crossed 
by the initial shock, there is still a significant amount of 
material attached to the core. Unlike S3, in S2 the obstacle 
mass-loss is faster initially, and when the obstacle is com- 
pletely shocked, the smaller amount of gas left in the outer 
layers, when detached, is not enough to increase the obsta- 
cle cross-section. To which extent a larger increase in the 
resolution could even change the results further in the same 
direction should be studied. 

The more efficient acceleration found in S3 and the 
higher values of the mass flux relies on the expansion of 
the obstacle gas after t ~ 4 x 10 5 s, as in SI. The expansion 
of the obstacle gas in S3 is very different from that of S2, 
as clearly seen comparing, e.g., the middle right panel in 



Fig. [8] and the bottom left one in Fig. [TT] When part of 
the expanded material in S3 enters the rarefied region sur- 
rounding the cylindrical tail behind the obstacle, it is accel- 
erated, reaching mean velocities v c ~ 10 9 cm/s. In S2, the 
expansion of the clump is not as important, so the pressure- 
driven acceleration has a smaller effect. 

4.3. Mass-loading 

On the axis, the acceleration of the obstacle material is 
basically linear with distance during the initial phase in 
all three simulations. A simple extrapolation of the linear 
growth results in 1.25 x 10 16 cm as the shortest distance to 
accelerate the obstacle gas to ~ Vj in SI. Taking the mean 
obstacle gas velocity, this distance becomes a factor of ten 
larger. In the case of S2 and S3, the extrapolation of the 
linear acceleration on the axis implies 6.5 x 10 15 cm, which 
becomes 3.25 x 10 16 cm for S2 and 1.3 x 10 16 cm for S3, 
taking mean obstacle gas velocities. However, we have seen 
that this process is critically dependent on the way in which 
the mass ablation is produced. If ablation mass-load is pro- 
duced at a smaller rate, the acceleration is more efficient 
for the ablated gas. In addition, if the extracted gas clumps 
widen enough to face the rarefied jet region, the obstacle 
material can be accelerated by pressure gradients up to a 
5-10% of the jet velocity already at ~ 10 14 cm, as for SI 
and S3. For S2, the obstacle gas does not expand enough to 
favor acceleration in the rarefied jet gas region. This ma- 
terial is confined around the symmetry axis, although this 
seems to be an effect of numerical dissipation. We remark 
that in either case, the dragged matter reaches z-velocities 
> v Q /x soon after being shocked, which confirms that it 
will remain in the jet as long as t& is at least few times 
shorter than t- } (see Section [2]). 

Efficient mixing of shocked jet and obstacle material 
could enhance acceleration (as in SI and S3), significantly 
reducing the distance needed by the entrained gas to reach 
the jet speed. When the obstacle gas expands due to heat- 
ing, efficient mixing may happen as a consequence of the 
growth of instabilities that disrupt the obstacle (SI), or 
due to the obstacle material expanding into a region oc- 
cupied by shocked jet flow (S3). However, if the mass- loss 
from the obstacle is too slow and/or the obstacle gas does 
not expand enough, all the ablated gas will flow through 
the tail. In this case, instabilities in the layer between the 
tail and the shocked jet flow would be necessary to favor 
mixing and subsequent acceleration of the obstacle gas. On 
the one hand, the tail is dense and surrounded by the rar- 
efied shocked jet gas, which is a very stable configuration 
in terms of Kelvin-Helmholtz or shear-layer instabilities. 
This is so because of the large relative inertia that the tail 
has under these conditions. On the other hand, any asym- 
metry in the medium surrounding the tail, which can be 
given by inhomogeneities in the jet or simply by the obsta- 
cle transverse motion, will result in the triggering of helical 
perturbations that are very disruptive. The latter cannot be 
observed in 2D axisymmetric simulations as those presented 
here, but should appear in 3D simulations. In conclusion, 
the distances given in the previous paragraph for obstacle 
material acceleration should be regarded as upper limits. 
The fact that they are of the order of the interaction z Q in 
this work, ~ 10 16 cm, likely implies that the jet will not get 
significantly diluted during the process and will be able to 
fully incorporate the obstacle matter. 
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Fig. 11. The upper panels show the same parameters as Fig. [2] for simulation S3 at time t^5x 10 4 s. The central and 
bottom panels are similar to the left ones in Fig. Hat t ~ 1.6 x 10 5 s, t ~ 4.3 x 10 5 s, t ~ 7.9 x 10 5 s, and t ~ 1.25 x 10 6 s 
from top to bottom and left to right. 
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Fig. 12. Axial cuts of different quantities for simulation S3 at t ~ 3.1 x 10 5 s, during the initial phase, before the 
shock completely crosses the obstacle. The left panel shows the normalized mass flux, normalized to its maximum value, 
3.7 x 10 22 g/s. The central and right panels show the mean velocity and density in the obstacle gas. 
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Fig. 13. Mass flux (left), mean density (central), and mean velocity (right) versus time in simulation S3 at z = 7.5 x 10 13 
(half grid, black lines) and 1.5 x 10 14 cm (end of the grid, red lines). 



The total mass flux in the simulated jets (considering 
the total jet cross-section) is M ~ 10 23 g/s ~ 10" 3 M o /yr. 
We have seen that, typically, the mass flux during the ini- 
tial phase is M ~ 8 x 10 20 — 3 x 10 21 g/s, which represents 
0.8-3% of the jet mass flux. However, during the develop- 
ment of instabilities and strong mass- loss, we found peaks 
of mass flux of M < 3xl0 22 g/s, i.e., 10-30% of the total jet 
flux. These episodes of dramatic mass-load could generate 
important deceleration and anisotropics in the jet flow (see 
also Sect. 15. Q . 

5. Discussion 

5.1. Mass-loading 

Jet formation likely involves a magnetocentrifugal mech- 



pairs as the jet matter content 


fe.fi:. iRuffini & Wilson 


19751: Blandfoi 


:d& ZnaiekLll977tl] 


Punslv & Coronitil 11990 


Beskin et al., 


19921 KomissarovL 


20041). Despite of being 



accretion-sustained, the magnetic field in the black-hole er- 
gosphere has a monopole-like structure that prevents ac- 
creted baryonic matter from entering into the jet. However, 
the jet may be already matter dominated a t some 10 3 i^sch 
from the black hole (e.g. lSikora et all 12005), probably wit h 



i P ly wit n 

an important content of protons (e.g. ICelotti et all , l2008f ) . 
In addition, the jet veloci ty, expected to increase in th e 
jet acceleration region (e.g. Komiss arov et all 120071 l2009f h 
may decrease at hig her values of the axial coordinate (e.g. 
iHoman et all l2QQ9h . All this can be explained by a jet ini- 
tially dominated by its Poynting flux, accelerated after- 
wards by magnetic- to-kinetic energy transfer, and finally 
slowed down due to mass-loading. 

Currently, the dominant mechanism and location of jet 
mass- loading are unknown. As shown in Sect. 14. 3| clouds 
and RG can effectively introduce baryonic matter into the 
jet innermost regions, with the latter acting as a sink (as 
long as it is not completely quenched by mass- loading; e.g. 
iHubbard fc Blackmanll2006D . A rough estimate of M can be 
done for the RG and BLR scenarios. In the former, adopting 
z Q ~ 10 17 cm, one may expect about o ne RG per year cross - 
ing the jet, delivering 10 28 - 10 29 g (jBarkov et all l2010h . 
i.e. M ~ 0.01-0.1 M cr (see Sect. [2]) for the jet considered in 
this work. In the BLR scenario, assuming a BLR filling: fac- 
tor ^ 10 -6 and p Q ~ 1.7 x 10 -14 g (see, e.g. lAraudo et all 
I2010D . one gets M - 0.1 M cr . For a leptonic lighter jet 
than considered here, with equal Lj but say Tj ~ 10, then 
M — >• M C r- Mass- loading can be important for more pow- 



erful jets if z Q is larger, since M ex z 2 (although this effect 
is partially compensated by the fact that M cr oc Lj). Note 
that in the case of stars, and neglecting the stellar wind- 
channel for mass- loading (i.e. RG or massive stars), z Q is 
limited by the tidal or the nuclear radiation intensity (see 
Sect. [I]). We conclude that, although a case-by-case analy- 
sis would be required, it seems clear that baryonic loading 
may be already important in AGN jets on pc scales. 



5.2. High-energy emission 

One obstacle covers a fraction of the jet cross section; 
about (R /Rj) 2 (= £) of the jet magnetic and kinetic en- 
ergy can be converted into internal energy, and a frac- 
tion of it, into non-thermal energy. The production of 
non-thermal accelerated particles in the considered re- 
gions can occur in shocks, turbulence, sheared flows or 
current sheets with reconnecting: magnetic fields (e.g . 



current sn e ets witn r econnecting magnetic nelas (e.g . 
ReesL 119781: iRieger fc Duffy! l2004HTammi fc DuffvL 120091 : 
Lvubarskv et al I l2Qlo( T These non-thermal particles, either 
protons or electrons, can interact with the ambient mag- 
netic, radiation and matter fields, producing high-energy 
emission. For protons, and in the context of jet-obstacle 
interactions, possible high-energy radiation channels are 
proton-proton (jpp) collisions in the entrained matter, pho- 
tomeson production if dense photon fields are present, 
and proton synch rotron for hi g h ma g netic fields and pro- 
ton energies ( e.g. Par fc LaorL Il997t lAraudo et all [2010; 
iBarkov et all I2010L l201ljT For electrons, relevant mecha- 
nisms would be synchrotron radiation at low energies, and 
synchrotron s elf- fSSC) and external Compt o n (EC) at high 
energ ies (e.g. iBednarek fc Protheroel 119971 : lAraudo et all 
120101 ). Secondary leptons, electrons and positrons (and also 
neutrinos) will result from pp interactions and photomeson 
production. At the jet base, the ambient photon fields can 
be also dense enough to absorb gamma rays, which will also 
result in e^-pair production. Therefore, processes of radia- 
tion attenuation and re processing cannot be neglected (e.g . 
iBlandford fc LevinsonLll995l : IBednarek fc Protheroel 1 19971 : 



Aha ronian et all 12008 ) , except for underluminous sources 
fe.g. IRieger et all 120081 : lAraudo et all . l2Q10t IBarkov et all 
120101 ). 

The simulations presented here show that for relatively 
homogeneous obstacles, the fraction £ can grow by more 
than one order of magnitude in few ta (~ 10 5 — 10 6 s), before 
the obstacle total destruction and mixing with the jet. It is 
seen in particular in Fig. in which about almost all the 
jet power within the grid is reprocessed, whereas the initial 
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Fig. 14. Spectral energy distribution beyond 1 GeV pro- 
duced by relativistic protons accelerated in the bow shock 
and interacting with the obstacle gas via pp. The accel- 
eration rate has been fixed to 0.01 qBc. The approximate 
sensitivities of Fermi and ground-based VHE instruments, 
for a source at ~ 20 Mpc, are also shown for comparison. 



£- value was only ~ 2.5 x 10 -3 (ho mogeneous grid). Th is was 
already predicted for instance in iBarkov et al.l ([2010h , and 
implies that the pp emissivity can grow very quickly while 
the obstacle pp cooling time is shorter than the dynami- 
cal time, assuming effective confinement of protons within 
the obstacle. When this cooling time becomes longer than 
the dynamical time, pp emissivity sharply decreases. This 
mechanism has been propos ed to explain t he da y-scale TeV 
flares observed in M87 (see lBarkov et al.L 120101 ). 

Figure [TH shows the expected pp gamma-ray spectrum 
for the times corresponding to Figs. [3] and i.e. the begin- 
ning and the £-peak phases, assuming that 10% of the pro- 
cessed jet power (the internal energy luminosity) in the grid 
goes to relativistic protons that interact with the shocked 
obstacle gas. As seen in the figure, the luminosity may grow 
by almost two orders of magnitude. In Fig. [I5j we show the 
evolution of the obstacle maximum radius enclosing den- 
sities > 10 10 cm -3 (at least in some regions), values for 
which the pp-cooling time is of the order of the dynamical 
time. For simplicity, we have taken a reference density of 
10 10 cm -3 , which yields a cooling time ~ 10 5 s. The fig- 
ures clearly show the difference in obstacle size evolution for 
the different cases discussed in the previous section, with 
a growth by a factor of six for SI, a factor of three for S3, 
and no growth for S2. For SI and to some e xtent S3, the 
result s are consistent with those presented in lBarkov et al.l 
(|2010f h which were obtained using an analytical approach 
for the obstacle dynamics. A more detailed account of the 
high-energy emission produced in RG-jet interactions will 
be presented elsewhere. 

Scenarios different from pp, like proton synchrotron 
emission, SSC and EC, have been also considered to ex- 
plain flaring as well as persistent emission from misaligned 
AGN and b lazars, in the context of jet - obstacle inter- 
actio n s (e.g. iBednarek fc Prothero"el Il997l : lAraudo et al.L 
120101 : IBarkov et all l201l[ ). In the case of proton syn- 
chrotron and EC, and also in photomeson production with 
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Fig. 15. Evolution of the obstacle maximum radius enclos- 
ing densities > 10 10 cm -3 (at least in some regions), for 
which the pp-cooling time is of the order of the dynamical 
time for SI (top), S2 (middle) and S3 (bottom). 



external fields, the decay phase of a flare would be longer 
than in the pp and SSC scenarios (as long as particle accel- 
eration runs). This is so because the former two processes 
are not so sensitive as the latter ones to the obstacle size, 
since the magnetic and external radiation fields may keep 
relatively constant during the event. 
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Interestingly, for power-law density profiles the obstacle 
ablated material is partially carried away in the form of a 
narrow tail behind a core, e.g. the RG itself. This, although 
to some extent an effect of the simulation resolution, indi- 
cates that in the RG case, the degree of detachment and 
homogeneity of the external layers of the star can affect the 
value of £, and thereby, the non- thermal luminosity. 

The quick obstacle acceleration by the jet can also 
have strong observational implications, separating two clear 
phases, one in which th e emission is rather isotropic, dur- 
ing th e first few td (e.g. lAraudo et al.L l2010b iBarkov et al.L 
l2010h . and one in which the obstacle remains move rela- 
tivistically but jet-driven bow shocks still surround them. 
In the latter situation, beaming effects would play a role, 
mak ing this phenomenon relevant also in the blazar context 
(e.g. IBarkov et all l201lh . 

Regarding a persistent component from interactions of 
many obstacles with the jet, its maximum luminosity can 
be estimated from N Q < (R Q /Rj) 2 > Lj < Lj, where 7V < 
(R /Rj) 2 > would be a typical jet cross s ection fraction 
cover ed by all the obstacles within z Q (e.g. lAraudo et al.L 
120101 ). We note that RG or occasional large clouds will be 
relatively infrequent, and thus contribute little to this quan- 
tity in the regions of interest (despite their occasional bright 
appearance). 

6. Summary 

Dense clumps and stars are expected to interact with/enter 
into the jets in the crowded and inhomogeneous central re- 
gions of AGN. Such events are likely to have strong con- 
sequences in the content, hydrodynamics and radiation of 
AGN jets. For typical jet luminosities and characteristics of 
the obstacle (BLR cloud, RG, etc.), significant amounts of 
mass are quickly dragged downstream, mixing with the jet 
flow. The mass injected into the jet can load it with baryons 
and reduce its Lorentz factor already on pc scales. The fast 
growth of the obstacles due to heating leads to a significant 
increase of the interaction cross section, enhancing the po- 
tential conversion of energy from the jet to non- thermal 
particles. This growth of the obstacle cross section strongly 
feeds back itself, and potentiates disruption. The accelera- 
tion of the obstacle by the jet, due to which the speed of the 
former approaches v } after several t&, changes the beaming 
pattern of the emission with time. Transient and persistent 
activity, related to different types of obstacles and their 
evolution within the jets, yields a rich and complex picture 
for the high energy emission of AGN, either misaligned or 
of the blazar-type. 
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